home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / zgbcon.z / zgbcon
Encoding:
Text File  |  2002-10-03  |  4.4 KB  |  133 lines
  1.  
  2.  
  3.  
  4. ZZZZGGGGBBBBCCCCOOOONNNN((((3333SSSS))))                                                          ZZZZGGGGBBBBCCCCOOOONNNN((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZGBCON - estimate the reciprocal of the condition number of a complex
  10.      general band matrix A, in either the 1-norm or the infinity-norm,
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE ZGBCON( NORM, N, KL, KU, AB, LDAB, IPIV, ANORM, RCOND, WORK,
  14.                         RWORK, INFO )
  15.  
  16.          CHARACTER      NORM
  17.  
  18.          INTEGER        INFO, KL, KU, LDAB, N
  19.  
  20.          DOUBLE         PRECISION ANORM, RCOND
  21.  
  22.          INTEGER        IPIV( * )
  23.  
  24.          DOUBLE         PRECISION RWORK( * )
  25.  
  26.          COMPLEX*16     AB( LDAB, * ), WORK( * )
  27.  
  28. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  29.      These routines are part of the SCSL Scientific Library and can be loaded
  30.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  31.      directs the linker to use the multi-processor version of the library.
  32.  
  33.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  34.      4 bytes (32 bits). Another version of SCSL is available in which integers
  35.      are 8 bytes (64 bits).  This version allows the user access to larger
  36.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  37.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  38.      only one of the two versions; 4-byte integer and 8-byte integer library
  39.      calls cannot be mixed.
  40.  
  41. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  42.      ZGBCON estimates the reciprocal of the condition number of a complex
  43.      general band matrix A, in either the 1-norm or the infinity-norm, using
  44.      the LU factorization computed by ZGBTRF.
  45.  
  46.      An estimate is obtained for norm(inv(A)), and the reciprocal of the
  47.      condition number is computed as
  48.         RCOND = 1 / ( norm(A) * norm(inv(A)) ).
  49.  
  50.  
  51. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  52.      NORM    (input) CHARACTER*1
  53.              Specifies whether the 1-norm condition number or the infinity-
  54.              norm condition number is required:
  55.              = '1' or 'O':  1-norm;
  56.              = 'I':         Infinity-norm.
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ZZZZGGGGBBBBCCCCOOOONNNN((((3333SSSS))))                                                          ZZZZGGGGBBBBCCCCOOOONNNN((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      N       (input) INTEGER
  75.              The order of the matrix A.  N >= 0.
  76.  
  77.      KL      (input) INTEGER
  78.              The number of subdiagonals within the band of A.  KL >= 0.
  79.  
  80.      KU      (input) INTEGER
  81.              The number of superdiagonals within the band of A.  KU >= 0.
  82.  
  83.      AB      (input) COMPLEX*16 array, dimension (LDAB,N)
  84.              Details of the LU factorization of the band matrix A, as computed
  85.              by ZGBTRF.  U is stored as an upper triangular band matrix with
  86.              KL+KU superdiagonals in rows 1 to KL+KU+1, and the multipliers
  87.              used during the factorization are stored in rows KL+KU+2 to
  88.              2*KL+KU+1.
  89.  
  90.      LDAB    (input) INTEGER
  91.              The leading dimension of the array AB.  LDAB >= 2*KL+KU+1.
  92.  
  93.      IPIV    (input) INTEGER array, dimension (N)
  94.              The pivot indices; for 1 <= i <= N, row i of the matrix was
  95.              interchanged with row IPIV(i).
  96.  
  97.      ANORM   (input) DOUBLE PRECISION
  98.              If NORM = '1' or 'O', the 1-norm of the original matrix A.  If
  99.              NORM = 'I', the infinity-norm of the original matrix A.
  100.  
  101.      RCOND   (output) DOUBLE PRECISION
  102.              The reciprocal of the condition number of the matrix A, computed
  103.              as RCOND = 1/(norm(A) * norm(inv(A))).
  104.  
  105.      WORK    (workspace) COMPLEX*16 array, dimension (2*N)
  106.  
  107.      RWORK   (workspace) DOUBLE PRECISION array, dimension (N)
  108.  
  109.      INFO    (output) INTEGER
  110.              = 0:  successful exit
  111.              < 0: if INFO = -i, the i-th argument had an illegal value
  112.  
  113. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  114.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  115.  
  116.      This man page is available only online.
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.